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Abstract 

The Goldstone theorem imphes the appearance of an ungapped mode whenever 
a continuous global symmetry is spontaneously broken. In general it does not say 
anything about the precise form of the dispersion relation nor does it imply that there 
is one massless mode for each broken symmetry generator. It is a well-established 
fact that even for relativistic field theories in the presence of a chemical potential 
Goldstone modes with quadratic dispersion relation, the type II Goldstone bosons, 
appear in the spectrum. We develop two holographic models that feature type II 
Goldstone modes as part of the quasinormal mode spectrum. The models are based 
on simple generalizations with U{2) symmetry of the well-studied holographic s-wave 
superfluid. Our results include Goldstone modes without broken generators but with 
unusual realization of symmetries and a frequency dependent conductivity of striking 
resemblance to the one of Graphene. 
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1 Introduction 



The AdS/CFT correspondence has proved to be a powerful tool to study strongly cou- 
pled quantum systems, with applications to QCD and condensed matter systems. Among 
the applications to condensed matter physics one of the most important examples is the 
construction of the so called holographic superfluids [U El 13] . 

The ingredients of such a model are a charged black hole describing a CFT at finite 
temperature and charge density. At sufficiently low temperature (or high chemical poten- 
tial) a charged scalar field develops an expectation value and triggers a symmetry breaking 
phase transition towards a superfiuid phase. In the simplest model the order parameter is a 
scalar and therefore we speak of an s-wave superfiuid. Furthermore it is possible to go to a 
decoupling limit in which the fiuctuations of the bulk metric are suppressed. The dynamics 
in this limit is completely specified by the fiuctuations of the gauge field and the scalar field 
[3]. Generalizations to p-wave [1] and d-wave E] superconductors, for which the order 
parameter has angular momentum / = 1 and / = 2 respectively are also known. Much of 
this development is reviewed in [71 [8] . 

The strength of the gauge/gravity correspondence is that it allows to study the real time 
dynamics of strongly coupled field theories rather easily. Linear response theory captures the 
behavior of a quantum system after an initial, small perturbation. It also applies to the late 
time behavior when an initially large perturbation has already sufficiently died out and enters 
the linear regime. The basic ingredient of linear response theory is the retarded Green's 
function. In the context of the AdS/CFT correspondence it was shown in [9l [10] how to 
calculate retarded Green's functions by imposing infalling boundary conditions at the black 
hole horizon. For black holes with non-degenerate horizons the retarded Green's functions are 
analytic in the upper half of the complexified frequency plane and have (an infinite series of) 
isolated poles in the lower half. These poles are the holographic quasinormal modes (QNM) of 
the black hole [HI US [131 [H] . Within the QNM spectrum, the ungapped modes play a special 
role, since they give the dominant contribution to the retarded Green's functions at low 
frequency and small momentum. Therefore they determine the hydrodynamic description of 
the system. In this paper we will study the quasinormal mode spectrum of a multi-component 
generalization of the simple, holographic, non-backreacted s-wave superfiuid with particular 
focus on the low lying and hydrodynamic modes. 

The study of the QNMs for the s-wave U{1) superfiuid was first carried out in [13]. 
Since the basic physics of superfiuids is the one of spontaneous symmetry breaking it can 
be expected that known results such as the existence of a Goldstone bosoiiQ carry over 
to the holographic models. Indeed, one of the main results of [15] was that the QNM 
spectrum in the superfiuid phase contains such an ungapped Goldstone mode with dispersion 
relation uj = ±Vsk + O^k"^). This mode can also be understood as the sound mode of 
the superfiuid and Vg is the sound velocitjj^ In the non-backreacted model these are the 

^It is also often called Nambu- Goldstone boson. For simplicity we will refer to it as Goldstone boson or 
Goldstone mode throughout the text. 

^ In [in] it was pointed out that this mode corresponds to the fourth sound. Whereas second sound 
is defined through temperature oscillations fourth sound is the phenomenon of sound propagation in the 
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only hydrodynamic modes in the broken phase. In the unbroken phase in contrast there 
exists a single hydrodynamic mode signaling the usual diffusive behavior of a normal fluid. 
Its dispersion relation is a; = —iDk"^, where D is the diffusion constant. The question 
what happens to this diffusive mode in the broken phase was also answered in [15]: it 
develops a purely imaginary gap u = —i'j — iDk"^. This is quite natural because the single 
purely imaginary mode can not move off the imaginary axi^ The hydrodynamics of the 
broken phase is fully captured by the Goldstone mode and the diffusion mode does the 
simplest thing it can to drop out of the hydrodynamic regime by growing the gap 7. Since 
this purely imaginary gapped mode has its origin in the universal diffusive mode of the 
unbroken phase we expect that it is a universal feature of a large class of superfluids, not 
only holographic ones. This mode will necessarily dominate the late time response in the 
order parameter to homogeneous perturbations and in regimes close but below the critical 
temperature where the gap 7 is rather small. Therefore the order parameter is bound to 
show a purely exponential decay towards its equilibrium value without any oscillation. In 
contrast for lower temperatures where 7 becomes large there are other low lying QNMs with 
real and imaginary parts in their frequency. In this low temperature regime the response 
in the order parameter is then an exponentially damped oscillation rather than a purely 
imaginary decay. This universal aspect of the late time response of superfluids was also 
emphasized in recent numerical studies of quenches of holographic superfluids in . 

In this paper we generalize the results on the QNM spectrum to models with U{2) 
symmetry. In the first model we simply add a second scalar field of the same mass, we will 
call this the ungauged model. A second model also includes gauge fields for the whole U{2) 
symmetry. The difference between the two models is as follows. In the ungauged model 
only the U{1) symmetry is local in the bulk. It has however a global SU{2) symmetry under 
which the scalar fields transform as a doublet. According to the holographic dictionary this 
model contains only one conserved current, corresponding to the single gauge field in the 
bulk. The dual field theory inherits of course the global SU{2) symmetry of the bulk but this 
symmetry is not generated by operators in the dual conformal field theory. This is similar to 
the decoupling limit in which we are working and in which the fluctuations of the metric are 
suppressed. The dual field theory has then strictly speaking no energy momentum tensor. 
In usual four dimensional Lagrangian field theories Noether's theorem guarantees that we 
can always construct a conserved charge generating a given symmetry of the Lagrangian. In 
holographically defined field theories the existence of a four dimensional Lagrangian is a priori 
not guaranteed and therefore Noether's theorem does not straightforwardly apply. This is the 
case here. Although the dual field theory has the SU{2) symmetry (and Poincare covariance) 
it does not contain operators generating these symmetries. We can speak of these symmetries 



superfluid component only 40J. This is the mode that survives the probe hmit in which propagation in the 
normal component of the dual fluid is prohibited. In this model there is no other sound mode so we will 
simplfy refer to it as the sound mode. Due to the nature of fourth sound it interpolates between second 
sound at T — Tc and normal sound at T = 0. 

■^Quasinormal modes are bound to come either in pairs aj„ and a)„ = '^^ fixed on the imaginary 
axis. This follows from rather generic symmetry considerations for retarded Green's functions, see appendix 
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as an outer automorphism of the operator algebra of the dual field theorjQ Physically the 
difference between the two models is that the ungauged one is a one-component fluid (there 
is only one notion of charge) whereas the gauged one is a two component fluid. In the latter 
case the charges are the expectation values of the zero- component of the currents in the 
Cartan subalgebra of the U{2) symmetry. 

Although this ungauged model does not contain conserved currents for the SU (2) sym- 
metry and therefore many of the standard proofs about existence of Goldstone bosons do not 
strictly apply we find a new ungapped mode in the QNM spectrum of the scalars. This mode 
is however not a standard Goldstone boson with linear dispersion relation but a so-called 
type II Goldstone mode whose energy depends quadratically on momentum. 

The second model we consider has a scalar field doublet coupled to the full set of U (2) 
gauge fields. We switch on a chemical potential only for the overall U{1) symmetry. Therefore 
the high temperature phase has the full U (2) symmetry. At low temperatures this symmetry 
is broken to U{1). In this model the dual field theory contains currents for all the U{2) 
symmetries. We can therefore also study the conductivities. 

In the context of condensed matter physics it has been pointed out long ago in [18] that 
such multicomponent superfluids have unusual Goldstone modes with quadratic dispersion 
relation. In the high energy context such models have been considered as models for Kaon 
condensation in the color-flavor locked phase of QCD in [T9l 120] again emphasizing the exis- 
tence of the quadratic Goldstone mode. Our gauged holographic model is a straightforward 
holographic analogue of the model in [121 120] and indeed we also find the presence of a 
Goldstone mode with quadratic dispersion relation. Let us also note that in the holographic 
context a type II Goldstone boson was found before in magnetized D3/D5 defect theory 

m- 

It seems useful to collect now some of the known theorems on Goldstone bosons (a very 
useful review on symmetry breaking and Goldstone modes is [22]). First we have of course 
the actual Goldstone theorem. Its proof assumes the existence of a conserved current 
such that the broken charge is Q = J df'xj^ (with d spatial dimensions). The theorem then 
states that spontaneous breaking of a continuous global symmetry implies the existence of 
a mode whose energy fulfills 

limw(fc) = 0. (1) 

The theorem by itself does not make any statement about the number of these modes, 
nor does it fix the ^-dependence of the frequency. In the presence of Poincare symmetry one 
can make however a stronger statement, namely that the dispersion relation of the Goldstone 
mode has to be linear and that the number of Goldstone bosons equals the number of broken 
generators. 

Lorentz symmetry might be absent however, either in principle such as in non-relativistic 
field theories or the system under consideration might be in a Lorentz symmetry breaking 
state, such as being at finite temperature or density. In these cases another theorem classifies 

string theory example for such a situation is the theory based on the small M — ^ superconformal 
algebra on the world sheet. This algebra possesses a large S0{4) — SU{2) x SU{2) symmetry acting on the 
four supercharges of which only one SU{2) is represented through chiral currents on the worldsheet. 
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Goldstone bosons as type I if their energy vanishes as an odd power of the momentum 
or as type II if their energy vanishes as an even power of the momentum in the zero 
momentum hmit. The number of type I and type II Goldstone bosons has to fulfill then 



rij + 2njj > Nbg , 



(2) 



where Nbg is the number of broken generators [23j. The number of type I and type II 
Goldstone bosons can be further constrained. Upon assuming that the broken symmetry 
generators obey {[Qa, Qb]) = Bab the number of Goldstone bosons has to fulfill [2ll EHl ET] 
(see also [22l EH ES] for more on counting rules of Goldstone bosons). 



We note that the gauged holographic model fulfills all these theorems. We have in total 
four symmetry generators. The symmetry is broken from U{2) to U{1) and so there are three 
broken generators. In the broken phase the charges corresponding to the overall U{1) and 
the Cartan U{1) generator inside SU{2) receive vacuum expectation values. Therefore the 
rank of the matrix B is two and so the number of type I and type II Goldstone bosons should 
add up to two. This is precisely what we find in the QNM spectrum, one ungapped mode 
with linear dispersion relation and one ungapped mode with quadratic dispersion relation. 

We also note that the ungauged model satisfies Goldstone's theorem and the counting 
rule of Chadha and Nielsen ([2]). It violates however the more refined counting rule (|3]). In 
a strict sense this model only has one symmetry generator since it has only one U{1) gauge 
field in the bulk. Therefore the counting rule ([s]) would suggest the existence of only one 
Goldstone boson, the number of broken generators is one and the matrix B vanishes trivially. 

This paper is organized as follows. In section [2] we briefly review a simple field theoretical 
model featuring type II Goldstone modes. This model has been introduced in the context of 
Kaon condensation in color-flavor locked QCD. It serves us as inspiration for constructing 
the holographic models. 

Section [3] is devoted to the analysis of the ungauged model. Since the well-known s-wave 
superconductor is a subsector of both the ungauged and the gauged model we also briefly 
review first the findings of [I5]. Then we show that even with this drastic simplification, i.e. 
not gauging the global SU{2) symmetry in the bulk, the model presents Goldstone modes 
with quadratic dispersion relation. Hence, within this model a type II NG boson is found as 
a consequence of having broken just one charge generator (the one associated to the U{1) 
symmetry) . 

In section |4] we study the fully gauged U{2) model. Then we analyze the fluctuation 
equations to linear order. They decompose into three decoupled sectors. One being the 
already encountered U{1) s-wave superfiuid, the other describing the non-Abelian sector in 
which the type II Goldstone mode resides and a third one with the unbroken U{1) symmetry. 
We proceed to study the conductivities which now arrange naturally into a two by two matrix. 
We show that the diagonal conductivities have delta-functions at zero frequency and are 
in this sense superconducting. Furthermore upon a change of basis we find a frequency 



ni + nil = Nbg - 




2 



(3) 
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dependence that is strikingly similar to the one of Graphene [22] • Furthermore we find 
indications that for temperatures below T = OATc another instability arises in the gauge field 
sector leading to an additional p-wave condensate. Then we study the low lying quasinormal 
modes and analyze the results. We find the type II Goldstone mode and also study the fate 
of the diffusion modes in the broken phase. Since now two symmetries participate there are 
two diffusion modes that in the broken phase pair up and can move away from the imaginary 
axis. We find that this is precisely what happens. Therefore the response in this sector does 
not show the purely exponential decay induced by the gapped pseudo diffusion mode of the 
[/(I) sector. 

We conclude in section|5]with a discussion of our results and an outlook to further possible 
studies of holographic type II Goldstone modes. 

Finally, in appendix [A] we present some general properties of matrix- valued Green's func- 
tions from which constraints on the quasinormal mode spectrum follow . In appendix |B] we 
collect technical details on how to actually compute the QNMs for coupled systems. 



2 A field theoretical model with type II Goldstone bo- 
son 

Motivated by the physics of Kaon condensation in the color-flavor locked phase of QCD 
the authors of [121 120] studied QCD at a nonzero chemical potential for strangeness. It 
was shown that at a critical value of the chemical potential equal to the Kaon mass, Kaon 
condensation occurs through a continuous phase transition. Moreover, a Goldstone boson 
with the non-relativistic dispersion relation w ~ appears in the Kaon condensed phase. 
To illustrate this fact, they considered the following (Euclidean) toy model: 

£ = (^0 + fi)<p\do - + d,<l)^di<p + M^<f)^<l) + A(0V)' , (4) 

where is a complex scalar doublet, 

/ ^. \ 

(5) 



As long as < M the masses of the four excitations in the model are the positive roots 
in u of 

{io ± = . (6) 

All are doubly degenerate. It is straightforward to check that at fi = M the global U{2) 
symmetry gets broken and the new vacuum can be chosen to be: 
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Studying the fluctuations of tlie doublet around this background one finds two massless 
and two massive modes with the following dispersion relations: 



= - 2M^ + 0(/) , (9) 

col =P^- , (10) 
ul=p'^ + 2fxuJi . (11) 

If we concentrate on the positive roots we see that cui is a normal, linear Goldstone mode. 



The positive root of equation (10) is 



c^3 = |^ + 0(p^). (12) 

This is the type 11 Goldstone mode. It has formally a non-relativistic dispersion relation. 
Since the underlying theory has however Lorentz invariance there is of course also a negative 
energy mode with quadratic dispersion. This arises as the negative root of 0J4,. Finally 0J2 
and are gapped modes with 

W4 = 2/i + 0(j92). (13) 

Since the symmetry breaking pattern is U(2) — )■ f/(l) there are three broken generators but 
only two massless Goldstone modes in the spectrum. This model fulfills all the counting 
theorems noted in the introduction. In particular the Chadha-Nielsen rule ^ is exactly 
saturated. The role of 074 is special. It is the mode that pairs up with the type II Goldstone 



mode in the dispersion relations (10) and (11). It has been argued that this mode is a 



universal feature and that its energy at zero momentum is exact and protected against 
quantum corrections |28l |30l |3T]. The spectrum obtained from this model is summed up 
in Figure [T] In our holographic models we will look for this special gapped partner mode 
of the type II Goldstone mode. It will turn out that the gauged and ungauged models 
differ significantly here: only the mode in the gauged model shows the characteristic linear 
dependence on the chemical potential. 

This simple Lagrangian model serves as our motivation and guideline to construct a 
holographic model featuring type II Goldstone modes. In fact we can use the same kind of 
matter Lagrangian in a holographic setup. According to the usual holographic dictionary 
a local bulk symmetry corresponds to a global symmetry in the boundary conformal field 
theory. We would therefore most naturally be led to a model in which we gauge the global 
U{2) symmetry of ^ and put it into an AdS Schwarzschild background. In order to trigger 
spontaneous symmetry breaking we introduce a chemical potential via a boundary value for 
the temporal component of the overall, Abelian U{1) gauge field. This is then our gauged 
model. 

Alternatively we might ask what are the minimal ingredients necessary to trigger sponta- 
neous symmetry breaking. The chemical potential resides entirely in the overall U{1) factor. 
The other three SU{2) gauge fields are not needed to achieve symmetry breaking. Therefore 
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0.5 1.0 1.5 2.0 M 



Figure 1: The spectrum of the field theoretical model. Below the critical value jj, = M 
there are four massive modes. The masses are M — fi and M + /i, the numbers indicate that 
they are doubly degenerate. In the broken phase ^ > M there are two Goldstone modes 
with exactly zero mass and two gapped modes. The special gapped mode has mass 2fi. 

we can choose as a sort of minimal setup a model in which the SU{2) symmetry stays global 
in the bulk of AdS. As already mentioned in the introduction this is a somewhat unusual 
realization of the symmetry from the boundary conformal field theory point of view. There 
are no conserved currents associated to this SU (2) symmetry, nevertheless all states and 
operators fall naturally into representations of this symmetry group since it is a global sym- 
metry of the bulk and it is also not broken by any of the boundary conditions. This setup 
constitutes our ungauged model and we will study it in detail in the next section. 

Let us note here one more technical detail: the field theoretic model of this section is most 
naturally viewed as living in four space time dimensions. In the following our holographic 
models will be dual to field theories living in three space time dimensions in order to stay 
as close as possible to the well-studied holographic U{1) s-wave superfiuid of [21 [T5]. This 
is however of no relevance to the essential features of the models, i.e. the existence and the 
nature of the hydrodynamic and Goldstone modes. 



3 The ungauged model 

We will now study the holographic model where the condensation of a charged scalar 
breaks a global SU{2) symmetry in the bulk. We shall look at the spectrum of quasinormal 
modes on both sides of the phase transition and study their dispersion relations. Since the 
simple U{1) s-wave holographic superfiuid constitutes a subsector of this as well as of the 
gauged model we will also use the opportunity to briefiy review the most salient features of 
its QNM spectrum. 

The minimal holographic model containing a type II Goldstone boson consists of a scalar 
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doublet of SU{2) charged under a f/(l) gauge field. The Lagrangian is given by 



(14) 



where 

*=(^), D, = d,-iA,, (15) 

and is the Abehan gauge field. The mass of the scalar field is chosen to be = —2/L'^. 
This is basically the same as the model in [2] except that we have added a second scalar field 
A with the same mass. Because of the degeneracy in the mass the model possesses in addition 
to the bulk-local U{1) symmetry a bulk-global SU{2) symmetry. Note that the global SU{2) 
symmetry is a priori not enough to set the field A(r) = 0. But we are interested in un-sourced 
static solutions for the scalar fields, i.e. we assume that the leading non-normalizable mode 
is not switched on. The solution space is then a two dimensional complex vector space 
spanned by the vevs of the operators dual to the scalar fields. On this parameter space we 
can act with the global SU{2) symmetry to set the operator corresponding to the field A 
equal to zero. Since now the non-normalizable and the normalizable mode of A are set to 
zero it follows that A(r) = 0. 

We will be working in the probe limit, in which the coupling of the gauge field is very large 
and the backreaction of the matter fields onto the metric can be neglected. The background 
metric is then taken to be the Schwarzschild-AdS black brane 

rfs^ = -f{r)dt^ + -f- + ^(rfx^ + d?/2) , 
/(r) 

M 

fir) = j-,--- (16) 

The horizon of the black hole is located at th = M^l'^l?!'^ and its Hawking temperature is 
T = In the following we use dimensionless coordinates 

( L^- \ 
(r, t, x,y) ^ thP, — t, — x, — y . (17) 

V th th J 

These rescalings allow us to set M = r// = 1 in the dimensionless system. In order to switch 
on a finite chemical potential in the boundary theory, the bulk Maxwell field 

A = x{p)dt, (18) 

must take a non-zero value at the boundary. 

The equations of motion for the background fields are 

X +-X ^X = 0, 19 

P f 

^"+(y + ^)^' + yJ^-yV' = 0. (20) 
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Notice that the system above is precisely the original U{1) holographic superconductor first 
studied in [2]. To ensure finiteness of the norm of the current at the horizon, we have to 
demand the scalar field to be regular whereas the gauge field has to vanish x{p = 1) = 0. 
With these boundary conditions, the asymptotic behavior of the fields at the conformal 
boundary is 

X = fi-- + o(^), (21) 



V-^ + ^ + o(l). (22) 

p p \py 

For the chosen value of the scalar mass, both terms in the scalar asymptotics correspond to 
normalizable modes [32]. Considering one or the other as the vacuum expectation value of a 
dual boundary operator leads to two different theories. In what follows we will consider only 
the case in which ipi is interpreted as the coupling and 1IJ2 as the vev of a mass dimension 
two operator. 

The dimensionless parameters are related with the physical quantities by 

a = ^1^, (23) 

" = ''''' 

where fi, n and Jo, (O) are the chemical potential, charge density and source and expectation 
value of an operator O of dimension 2, respectively. From now on we set L = 1. In the 
following we will work in the grand canonical ensemble. In practice we vary the dimensionless 
parameter /2. Because of the underlying conformal symmetry this can then be thought 
of as either fixing the chemical potential /i and varying the temperature T or fixing the 
temperature and varying the chemical potential. We define the temperature by T/Tc = Pc/p 
and fix /X = 1. 

Spontaneous symmetry breaking is driven by low temperature or high chemical potential. 
It triggers a non trivial expectation value for the scalar field without switching on any source 
Jo- For small p, the scalar field is trivial and the gauge equation is solved by x = p{^ — 1/p) 
and ip = 0. The system is then in the symmetric phase. However, by decreasing the 
temperature the system becomes unstable towards condensation of the scalar [U [2] . In [15] 
it was shown that at the critical temperature indeed the lowest quasinormal mode of the 
scalar field becomes unstable, i. e. it crosses over to the upper half plane. 

The free energy density of the system is given by the on-shell renormalized action, 

F = -TSren = ^/in - dr . (27) 
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The second term vanishes in the absence of a condensate and it works against the phase 
transition if it is present. In Figure |2] the free energies for the symmetric and broken phase 
are compared. It is clear that for T < Tc the condensate solution is always preferred and 
therefore the system undergoes a second order phase transition to the superconducting phase. 
Note that the presence of the second scalar plays no role for the phase structure. It simply 
vanishes in the broken and unbroken phase A = 0. 



< o, >^ 

F 

T* 




Figure 2: (Left) The free energy of the trivial (blue) and condensate (red) background so- 
lutions at low temperatures, T < Tc. (Right) Value of the condensate in the grand canonical 
ensemble as a function of T/Tc. 



In order to extract the quasinormal mode spectrum, we switch on fluctuations of the 
background fields 

= i7]{p,t,x),ijip) + a{p,t,x)) , (28) 
^ = ixip) + (^t{p,t,x))dt + ax{p,t,x)dx . (29) 

We do not include transverse fluctuations because they decouple from the interesting physical 
features of the model at hand. 

In the normal phase, i.e. expanding around ipip) = 0; the system reduces to the U{1) 
holographic superconductor studied in [TFJ with two copies of the scalar fluctuations, 

/." + .'(/'+f) + (<^-^-m^).^0. (30) 

„ 2f , k'^ ook 

fa^ H -at -a^ =0 , (31) 

p p'' 

2 7 

K + /'< + y«. + y«t=0, (32) 

iuj , ik , , , 

ya; + -<=0, (33) 
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where s stands for both t] and a fluctuations. The equation for the complex conjugate 
scalar s can be obtained by changing the sign of the potential x (30). The frequency and 
momentum are related to the physical ones by 
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(34) 



CO = 



AttT 




(35) 



The scalar and gauge fluctuations completely decouple in the symmetric phase. This is 
a consequence of working in the probe limit. The quasinormal mode spectrum of the U{1) 
field in the normal phase is just that of an electromagnetic field on an AdS-Sch background. 
The longitudinal fluctuations contain one hydrodynamic mode, u = —iDk"^, reflecting the 
diffusive behavior of normal fluids. In physical units D = 3/{AttT). Due to the lack of an 
energy- momentum tensor for the dual field theory in the probe limit, the diffusion pole is 
the only hydrodynamic mode in the unbroken phase. 

There are two copies of the scalar fluctuations. The quasinormal modes of t] and a move 
towards the origin when decreasing the temperature, whereas the modes of fj and a have 
larger masses and widths the smaller the temperature. As we approach the critical temper- 
ature T = Tc, the lowest quasinormal modes of rj and of cr become massless, triggering the 
phase transition: the scalar field acquires a non trivial vev in order to avoid its fluctuations 
to become tachyonic. By symmetry we can choose the condensate to reside completely in 
the ijj field. The fluctuations a couple then to the gauge field fluctuations just as in [15] . 
Therefore the QNM spectrum in this sector contains a Goldstone mode with linear dispersion 
relation u = ±Vsk + 0(A;^). This is the usual type I Goldstone boson associated with the 
breaking of the gauge U{1) symmetry. As shown in [T5] it can be interpreted as the sound 
mode of the dual superfluid in the broken phase. What happens then to the QNMs in the 
fluctuations of the second scalar 77? At the critical temperature there is also an ungapped 
mode present since its QNM spectrum is simply another copy of the scalar sector. Since 
there are no operators generating the SU{2) symmetry in the dual field theory standard 
arguments about the appearance of Goldstone modes do a priori not apply. Three logical 
possibilities arise then: the mode could become unstable for T < T^, it could become gapped 
again or it stays ungapped, playing the role of an unexpected Goldstone boson for the broken 
bulk-global SU{2) symmetry. Shortly we will see that the last possibility is realized and that 
the massless mode of 77 will indeed correspond to a type II Goldstone boson with quadratic 
dispersion relation, u oc k'^. 
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In the broken phase, the equations of motion read 



o = fv" + v'{f' + j] + 

2/ 



= f6" + 6'if' + 



= /C" + C' \^f' + 



o = K + 
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-at 
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+ 2ip ] at -a^ - 2iuip6 - ixi^C ■ 



= /< + /'< + (y - 2V^) «x + y at + 2iki;5 . 
0= ^a; + ^< + 2^'5-2#', 



(36) 

(37) 
(38) 

(39) 

(40) 

(41) 



where we have divided a = ( + i5 into real and imaginary part. The system (37)-(41) is 
again the one studied in [T5j . This sector, that also appears in the gauged model that will be 
presented afterwards, decouples from the additional scalar fluctuation rj. Notice that even 



if (36) is formally the same as in the normal phase, the background x is different leading to 



non trivial features in the rj sector such as the presence of a massless excitation. 





Figure 3: Real (left) and imaginary (right) parts of the lowest scalar QNMs as a function 
of the chemical potential. Solid lines correspond to the unbroken phase. For the broken 
phase dashed lines stand for modes of the additional scalar while dotdashed lines represent 
the modes common to the U{1) holographic superconductor. 

Figure [3] shows the spectrum of quasinormal excitations of the scalar doublet. In the 
normal phase we have two degenerate copies of the spectrum that partially split after the 
phase transition. It is clear that the two lowest excitations become massless at the critical 
chemical potential and then remain massless in the superconducting phase. They can be 
identified with the two Goldstone bosons at the phase transition. The rest of the excitations 
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remain gapped in the broken phase. Notice that the first f] excitation (dashed black hne in 
figure [3]) does not follow the expected universal behavior in the broken phase, i.e. it is not 
linear in fi. This mode is the equivalent of the special gapped mode U4 in the field theoretical 
model of section |2} However, it has already been mentioned that the ungauged model does 
not satisfy all the theorems about symmetry breaking and therefore deviations from the 
universal behavior should not be surprising. The behavior of the gapped modes is actually 
similar to that of the U{1) model modes. In the unbroken phase we can distinguish the 
modes that come from the s-type of fiuctuations from the ones that come from the complex 
conjugate s fiuctuations. The former become lighter whereas the latter become heavieiQ In 
the broken phase it is more useful to use real and imaginary parts, at least for the scalar 
that mixes with the gauge fields fiuctuations, i.e. the lower component of the scalar in our 
conventions. So we can not a priori talk of s and s type fiuctuations. We still can study 
to which modes the s and s type modes connect to in the broken phase. Here we see an 
interesting pattern: the s type modes split in the broken phase whereas the s type modes 
stay almost degenerate close to the phase transition (at least at zero momentum). This is 
surprising given the fact that the fiuctuations correspond to two completely different systems, 
one coming from a single differential equation whereas the others come from a complicated 
system of coupled equations. However, for small temperatures they split and actually the 
real part of the lowest one for the U{1) sector goes to zero at a finite temperature. For 
temperatures below T 0.63 Tc it becomes a purely imaginary mode. 



Sound mode: There are two massless modes in the broken phase. The first one is the 
type I Goldstone boson appearing because of the spontaneous breaking of the U{1) gauge 
symmetry. In [15], it was shown that this mode corresponds to the sound mode of the dual 
superfiuid and that in the hydrodynamic limit it has a linear dispersion relation 

uji = ±{vsk + h k^) - iV, k"^ , (42) 

where Vg is the speed of sound and Vg is its attenuation. It turns out that the quadratic 
part of the dispersion relation also has a real component. This component is very small and 
subleading compared to the linear term that determines the speed of sound. In [15] this real 
quadratic part has not been studied. 

For very small temperatures the velocity approaches its conformal value v1 = 1/2 while 
the width goes to zero, see figure |4j Close to the phase transition, the speed of sound has a 
mean field behavior as a function of temperature 

v'g ^ 2.8 (1 - I) • (43) 

As expected, at T = the speed of sound vanishes. This can be traced back to the fact 
that at the phase transition the total mass ml = — /x^ fulfills ml = f ^ = 0, as expected, 
and hence the complex scalar field, charged under a U{1) symmetry, becomes massless. 

'^This behavior is reversed if we had taken the chemical potential to be negative. 
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Figure 4: Speed of sound and damping for the sound mode. The speed of sound goes to 
zero at the critical temperature. The damping constant first rises quickly and then falls off 
again. Precisely at the critical temperature its value is such that the sound modes connect 
continuously to the scalar modes that become massless there. The peak in the damping 
constant sits close to the critical temperature and was not resolved in [T3] . 



Indeed, one can write down the effective action, analogous to Q, for a complex scalar 
field with mass M, in the presence of a chemical potential for a U{1) symmetry that is 
spontaneously broken. The excitations on top of the [/(l)-breaking background have a 
dispersion relation equal to (|8])-([9]), being (|8| the type I Goldstone boson. It is a general 
feature of these linear sigma models that the coefficient in front of the linear term in the 
momentum depends on ml, as can be explicitly checked for the case at hand (see ([8|). 
Therefore, at the phase transition the leading term in the dispersion relation is of 0{k'^); 
this effect can be seen very clearly with numerical methods, as shown in Figure [5] Since 
the quasinormal mode spectrum has to vary continuously through the second order phase 
transition the real and complex coefficients of the fc^ term have to coincide at T = with 
the ones obtained from the massless scalars in the unbroken phase. Numerically we find 
b{Tc) = 0.22 and r,(Tc) = 0.071. 



Pseudo diffusion mode: In the unbroken phase our model has only one hydrodynamic 
mode, the diffusion mode u = —iDk^ + 0{k'^) with D = 3/(47rT) in physical units. The 
shear and normal sound modes have their origin in the metric fluctuations and therefore are 
absent in the decoupling limit we are studying. The phase transition to the broken phase 
is second order. For the spectrum of quasinormal modes this implies that the modes of the 
broken and unbroken phase must connect continuously through the phase transition. In 
the case of the diffusion mode there must therefore exist a quasinormal mode with purely 
imaginary frequency. Hydrodynamics implies however that the only ungapped modes are 
the sound modes corresponding to the type I Goldstone mode. Not too far from the phase 
transition, i.e. for T < Tc the diffusion mode of the broken phase must develop into a mode 
with dispersion relation 

uj = -i-f{T) - iD{T)e + 0{k^) , (44) 
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Figure 5: Dispersion relations of Re a; (left) and Imo; (right) at T = Tc for the type I 
Goldstone boson in the system studied by [15]. The behavior Rew ~ becomes quadratic 
right at this temperature: Recu = bk"^. The coefficient is 6 = 0.22, which in turn matches 
the value that one finds if approaches Tc from above (i.e. from the unbroken phase). 



as shown in Figure [6j 

We might say that the diffusion mode develops a gap in the broken phase and becomes 
what has been called a pseudo diffusion mode in [TJj. Precisely at zero momentum k = 
this gapped pseudo diffusion mode explains the observation made in [17] on the late time 
response of holographic superconductors. For temperatures T < Tc the pseudo diffusion 
mode is the mode that lies closest to the real axis and therefore it dominates the long time 
response to any perturbation, such as the quenches studied in [17]. It follows that the order 
parameter shows a purely exponential decay since this mode does not have a real frequency. 
The existence of that mode can ultimately be traced back to the universality of the diffusion 
mode in the unbroken phase. We expect therefore the pseudo diffusion mode to be a universal 
feature of a wide class of superfluids (not necessarily holographic ones). 

The gap 7 grows as the temperature decreases. On the other hand there are quasinormal 
modes (connecting to the QNMs in the scalar sector of the unbroken phase) whose imaginary 
part is only weakly dependent on the temperature. At a certain crossover temperature T, 
the gap of the pseudo diffusion mode is bigger than the imaginary part of these modes, as 
shown in Figure [7] Then the response pattern changes from a purely exponential decay to 
an exponentially damped oscillation. Numerically we find that the crossover temperature is 

= 0.69 Tcj^ Such crossover changes in the long term response appear frequently in the 
details of the quasinormal mode spectrum of holographic field theories, [331 EH |35] . In fact 
this purely exponential decay applies not only to the order parameter but to all operators 



that correspond to the fields participating in the fiuctuation system (37)-(41), e.g. charge 
density or x-component of the current. 

For finite momentum the response pattern is expected to be different however. Now one 
also has to take into account the sound mode. While precisely at zero momentum the sound 



^This is lower than in the model of [17]. The difference is presumably due to the fact that we work in 
the decoupling limit. 
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Figure 6: (Left) Dispersion relation of the gapped pseudo diffusion mode in tlie broken phase 
for three different temperatures. The gap widens as the temperature is lowered. (Right) Gap 
7 as a function of the reduced temperature T/Tc. As one approaches the critical temperature 
from below the gap vanishes linearly. 




Figure 7: (Left) Continuation of the second and third scalar QNM into the broken phase. 
The real part grows as the temperature is lowered whereas the imaginary part shows very 
little dependence on T. (Right) The gap 7 (blue line) and the imaginary part of the lowest 
(scalar) mode fluctuation (red line) in the broken phase are shown as function of T/T^. 
At ^ 0.69 Tc the imaginary parts cross. For lower temperatures the late time response 
is not dominated anymore by the pseudo diffusion mode and consequently is in form of a 
exponentially decaying oscillation. 
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mode, i.e. the Goldstone mode, degenerates to a constant phase change of the condensate at 
small but non-zero momentum the long time response should be dominated by the complex 
frequencies (42). If one looks however only to the response in the gauge invariant order 
parameter \0\ the Goldstone modes, being local phase rotations of the order parameter, are 
projected out. 



Type II Goldstone mode: The second massless mode is the Goldstone boson associated 
with the breaking of the bulk-global SU (2) symmetry. It can be fit to a quadratic dispersion 
relation of the form 

uii = ±bk^- ic + 0{k^) , (45) 

in the long wavelength limit. Therefore it has the characteristic of a type II Goldstone mode. 
In Figure [8] the dispersion relation for the t] massless mode is shown for various temperatures 
as well as its fit to the hydrodynamic form. It is clear that there is a good agreement in the 
regime of validity of the low energy limit. 



Re toil 



Im tt)l; 





Figure 8: Real (left) and imaginary (right) parts of the type II Goldstone mode as a 
function of the momentum for T/T^ = 0.9998 (blue) and T/T^ = 0.704 (red) . The solid 
lines correspond to the numerical result while the dashed lines are the quadratic fit to a 
dispersion relation uu = hk'^ — ick"^. 



The coefficients in the hydrodynamic dispersion relation (45) as a function of the tem- 



perature are shown in Figure |9| Close to the phase transition they have a linear dependence 
in the reduced temperature 



h{T) 



c{T) 



0.22 + 0.049 1 - 



(46) 



0.071 - 0.0014 1 - 



T 



near 



Notice that at the phase transition the sound mode and the type II Goldstone must behave 
in the same way due to continuity of the modes through the phase transition and the fact 
that they are degenerate in the normal phase. In fact, at the transition b = b = 0.22 and 
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c = = 0.071, values that of course coincide with those of the lowest scalar mode in the 
normal phase. On the other hand, it is interesting to notice that in the broken phase the 
behavior of the coefficients of the type II Goldstone is completely different from that of the 
coefficients of the sound of the superfiuid. Unlike the sound velocity, that vanishes at the 
phase transition, the coefficient b of the type II Goldstone mode takes a finite value at the 
critical temperature. This result of course persists for the gauged model. The attenuation on 
the other hand, as it happens for the U{1) sector, has a finite value at the phase transition 
and then decreases with temperature, reflecting the fact that the fluid is more ideal the lower 
the temperature. 



b 




0.4 0.5 0.6 0.7 0.8 0.9 1.0 T, 



Figure 9: Coefficients of the type II Goldstone mode dispersion relation oju = bk"^ — ick"^, 
as a function of the temperature. Dependence with temperature is very mild. 



4 The gauged model 

Let us now discuss the fully gauged model. Consider the following Lagrangian for a 
complex scalar field living in the fundamental representation of a f/ (2) gauge symmetry of 
the bulk. 



where 




and c = 0, 1, 2, 3 is the color index. The field \1/ plays the role of the condensate. The 
expectation value of its dual operator thus triggers the spontaneous breaking of the U{2) 
global symmetry of the boundary theory. For simplicity, we set A = in the background. Tc 
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are the generators of U{2): 



{Ti,Tj} = ^6ijl, 



-I, 



T, 



2 
1 



{To,T,} = -ai 



(49) 



Notice that we are again working in the probe hmit, so the background metric is taken 

the ol 

e(r) 



to be the Schwarzschild-AdS black brane of (16). On the other hand, the gauge field is now 

I (3) 



A, 



(0) 



$fr) 



(50) 



The rest of the components of the gauge field being zero. As in the previous section, we will 



use dimensionless coordinates defined by the rescaling given in (17) 
The equations of motion for our ansatz are 



r + 



f P 



4/2 



m 

T 



^2 



+ ^($-6) =0, 

P 2/ 

e" + -0' + ^(<i>-0) = o. 

P 2/ 



(51) 
(52) 
(53) 



Notice that from (53) it follows that we can not simply switch on $ without also allowing for 



a non-trivial G. We are of course only interested in switching on a chemical potential in the 
overall U{1), and therefore we will impose 0(p — )■ oo) = and allow for a finite boundary 
value of 

The coupled system of equations above can be simplified by defining x = | ($ — ©) and 
^ = ^ ($ + 6). Using (52) and (53), we see that the resulting equations for these fields ar^ 



f 



2 2 

P f 

„ 2 , 2^2 



X 



P 



-X 



f 



-X 



e + k' 
p 



0, 
0, 

: 



(54) 

(55) 
(56) 



where we have redefined tp — )■ ^/2'^. As usual we choose the boundary conditions x(p = 1) = 
0, ^(p = 1) = along with regularity of \E'. Having a dual field theory with only one finite 
chemical potential switched on, implies that x and ^ must take the same non trivial value 
at the boundary in order to ensure that G vanishes asymptotically. Notice that ^ decouples 

^These equations of motion correspond to the probe limit of the system studied in j36j as a dual of 
superconductors with chemical potential imbalance. Notice however that in [35] the gauge symmetry was 
U{1) X U{1) instead of U{2) as in the present setup. 
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completely. The remaining system (54)- (55) is again the background found for the widely 
studied s-wave U{1) holographic superconductor. Therefore, the background of the U{2) 
gauge model contains the Abelian superconductor plus a decoupled conserved U{1) sector. 

The field x li^s in the direction of one of the broken generators, which is the linear 
combination ^{T^ — Tq), whereas ^ lies in the direction of the preserved U{1) given by 
K^3 + To). 

The asymptotic expansion of the fields near the conformal boundary reads 



X = Px 



- — + 

p 

P 



P P^ 



(57) 
(58) 
(59) 



The map of the various coefficients in the previous equations to the boundary conditions is 
fi^ = Ji^ = Ji. We will again focus in the O2 theory exclusively, henceforth we will demand 
^1 = 0. 

Equations (54)-(55) allow for solutions with a non- vanishing condensate, and therefore 



\{T^ — To) will be spontaneously broken. This solution must be found numerically, since the 



system is non-linear. However, (56) does have an analytic solution 



P 



(60) 



and thus n^ = p. 

When the symmetry is not broken, = 0, the equation for x has of course 



X = ^ 1 



as a solution as well. Therefore, in the unbroken phase 



6 = 0, 



$ = 2/i 1 



P 



(61) 

(62) 
(63) 



This behavior reflects the fact that T3 is completely independent from Tq in the unbroken 
phase. However, once we switch on the condensate, the interplay between T3 and Tq (recall 
that the remaining symmetry is a combination of the two) makes it impossible to set only 
one of the fields to zero. 

Finally, let us mention that in order to relate the dimensionless parameters with the 



physical ones, we need to apply the same dictionary (23)-(26) used for the ungauged model. 
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4.1 Charge Density in the broken phase 



According to [HI |2l] we can expect the presence of type II Goldstone modes if the broken 
symmetry generators fulfill 



(64) 



with at least one Bah 7^ 0. In our case we have [Qi,Q2] = iQz- Therefore in the broken 
phase we are interested in a non-vanishing expectation value for the charge density operator 
(Qs) = n^. As we argued previously, in the unbroken phase we necessarily have 9(r) = 0. 
This happened since both x ^"^^ ^ obey the same differential equation and the integration 
constants had to be set equal in order to do not switch on a source for O. Now we would 
like to find out whether or not an expectation value for O will be spontaneously generated 
in the broken phase. 

Independently of the phase the field associated to the unbroken combination of generators 
is given by (60). Since 6 = ,^ — then 



(65) 



Hence, what we want to check is the difference between the leading and the subleading 
coefficients of x as a function of the temperature. The numerical result is shown in Figure 

m 



Figure 10: Charge density of 9, uq, as a function of the temperature T/Tc. 



So we conclude that precisely at T < Tc this difference is switched on and an expectation 
value for (Qs) appears. This can be taken as a clear indication for the appearance of type 
II Goldstone bosons in the spectrum. 
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4.2 Fluctuations of the gauged model 



In order to study the quasinormal spectrum and the conductivities of the system, we 
switch on longitudinal perturbations on top of the background, so that 



= {r]{t,p,x),<il{p) + a{t,p,x)), (66) 

= mp) + a'f\t,p,x))dt + a'^^\t,p,x)dx, (67) 

A(^) = a^t^\t,p,x)dt + a^^\t,p,x)dx, (68) 

^(2) = aP{t,p,x)dt + a'^^\t,p,x)dx, (69) 

= {Q{p) + a't\t,p,x))dt + a'^\t,p,x)dx. (70) 



4.2.1 Perturbations in the Unbroken Phase 

In the normal phase, the background value of the condensate vanishes. Moreover, we 
have Q{p) = 0. The equations of motion for the perturbations read 

+ + + = 0, (71) 



/ P 











fp' 




I 2 /(c) 

+ -at - 
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a(^) ■ 
fp' ^ 




f 




uk ( 
+ ^a, 




^ /(c) 

T' 


+ ^4 
P^ 



0, (72) 
0, (73) 
0, (74) 

where s G {f], a}. Since the color indices do not see each other the system is the same one as 



(30)-(33) except that there are four copies of the gauge field fluctuations. Due to the chosen 
normalization of the U{2) generators the gauge field background $ enters with an additional 
factor I compared to (30). The quasinormal mode spectrum is the same as the one of the 



holographic s-wave superconductor except that the scalar modes are doubly degenerate 
and the gauge field modes are fourfold degenerate. In particular there are four copies of the 
hydrodynamic diffusion mode u = —iDk"^. 



4.2.2 Perturbations in the Broken Phase 

The equations of motion in the broken phase decouple in two sets: one mixing the (0) — (3) 
colors of the gauge field and a fluctuations and the other mixing the (1) — (2) colors and the 
r] fluctuations. 

Writing a = C, + i5, the equations of the (0) — (3) sector are 
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(3) _ (0) 

' ^^^^ 

= /af °) + ^af ) - f vl.^ + - ^af + v&^af ^ - 4C^X - 2^u^5 , (77) 

P V P J P 

= K(°) + /'af ) + - v[/2) 4°) + ^af) + v^M^) + 2. W , (78) 

= /af ) + ^af ) - f + ^) af ) - ^af + vl^^a^^ + K^X + 2.c.^5 , (79) 
P \ P^ J P^ 

= fal^^) + faf^+{^-'^^^a^^^ + '^af'^ + ^^af-2xk5^, (80) 

= !|af ) + ^af ) + 2vl/'5 - 2VI/5' , (81) 

= !|af ) + !^af - 2M/'5 + 2M/y . (82) 



It is trivial to show that by defining new fields = \{a!i^ ± af^) and aii^^ = ^(ai^'* ± 
the system further decouples into a coupled system for the scalar fiuctuations and 
a^T'' and a background independent set of equations for the U{1) gauge field o!i^\ The 



first subsystem reproduces the eoms (37)-(41) and therefore corresponds to the s-wave U{1) 
superconductor contained in the U{2) model. On the other hand, the field a^"^^ corresponds to 
the preserved gauge symmetry surviving the U{2) — )■ U{1) spontaneous symmetry breaking. 
The quasinormal mode spectrum in this sector is therefore the same one as in [15] plus the 
QNMs that are stem from a f/(l) gauge field in AdS^^. In particular the hydrodynamic modes 
in this sector are the sound mode and the diffusion mode of the unbroken f/(l). 

From now on we will concentrate on the remaining fields. We will call this remaining, 
inherently non-Abelian sector the (1) — (2) sector and will show that the expected type II 
Goldstone boson resides there. Writing r] = a + i/S, we find the following equations in the 
(1) - (2) sector: 
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fa + + j« + _ ;3_ 

( + V (83) 



2p2 - ' 2/ * ; 2/ 



= /af) + ^af)- fvI/2 + ^^a«-^a«+^e44^)-2$vl>a-2zc^v^^ (85) 
P \ P^J P^ P^ 



+2ifc^^ , (86) 



o I t 9 a:; ^ 9 

P VP/ P P 



= /af) + /Vp)+(^-^^ + ^)ai^)+2z^a«+.e^ain^^^^^ 



+2i/c\l>a , 



= ^af) + ^af) + i(af e-afe')+2vl/'/3-2/3'vl>, (89) 







^«^(2) + J^'m _ 1 (^^;(i)0 _ ^DQ'j + 2vl;'a - 2a'vl; . (90) 



A comment is in order here. This system of equations could be written in a more compact 
form by using complex field variables rj and a[^l ± ia[^]. . One has to keep in mind then that 
the field equations one needs to solve for the QNM spectrum for the complex conjugate 
fields are not the complex conjugate equations since one has to demand infalling boundary 
conditions on the fields and on the complex conjugate fields simultaneously. This aspect is 
somewhat clearer if one works with the (formally) real field variables on paying the price of 
writing a somewhat lengthy system of equations. 

Up to linear order in perturbations, there are three decoupled sectors in the system. 
Two of them belong to the '(0) — (3) sector' and they are a copy of the U{1) holographic 
superconductor, already extensively studied, and the preserved U{1) gauge symmetry. The 
main features of the spectrum of this sector have already been presented in section |3] since 
it is also a subsector of the ungauged model. On the other hand, the so called '(1) — (2) 
sector' has not been studied before. The physics in this sector is quite different from the 
holographic superconductors studied up to now and we will concentrate on it in the rest of 
this paper. 
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Before studying the quasinormal modes we will focus on a simpler problem, namely the 
conductivities. 



4.3 Conductivities 

In order to study the conductivities via Kubo formulae, it is enough to solve the linearized 
equations in the limit k = 0. The retarded correlators that we are interested in have the 

form Qn ~ {^<^fc)' "^(c')) ' "^i^h c, c' color indices. 

We will be applying the prescription of pi] for computing Green functions in the presence 
of holographic operator mixing. If one has a set of fields $/, the two-point correlation 
functions will be 

gij = hm (AM^f j(A)' + Bjj) , (91) 

A— >-oo ^ 

where the matrix J-fc(r) is nothing but the bulk-to-boundary propagator for the fields, nor- 
malized to be the unit matrix at the boundary. The matrices A and B can be read off from 
the on-shell renormalized action. The corresponding DC conductivities are given by the 
following Kubo formula 

ajj=hm( -Quito, 0)] . (92) 



At vanishing momentum the longitudinal components of the gauge field perturbations 
decouple from the scalar perturbations, as well as from the temporal components of the 



gauge fields. Moreover, the constraints (eqs. (81 82) and(89 90)) become trivial. Since we 



know that the system splits into the (0) — (3) and the (1) — (2) sectors we can rearrange the 
a^x^ fields in two vectors 

€(0~3)(P) = («f (P),«f (P)) and $r(i_2)(p) = {a^x\p),a^x\p)) ■ (93) 
One can check that in our case the A, B matrices take the simple form 

fir) 

A=-^-^l, B = 0, (94) 

for both sectors. A priori we would have a 4 x 4 matrix of conductivities. We know however 
that the fluctuations in the (0) — (3) and the (1) — (2) sector decouple from each other. 
Therefore we can restrict ourselves to study two independent 2x2 matrices of conductivities. 

4.4 Conductivities in the (0) — (3) sector 

The = equations of motion for ai^"* and aH^ can be simplified by using the already 
defined a!x^ and ai^-* fields. This results in 

^ K-./a-.fo^), (95) 
= K« + /Vj-) + (y - 2vI/2^ 4-) . (96) 
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We see that the resulting system of equations is now completely decoupled. We only have 

two diagonal conductivities cr++ and a , corresponding to the unbroken f/(l) diffusive 

sector and a mode which is associated to the broken f/(l) coupling to the condensate. The 
former is the same as in the unbroken phase and of no further interest for us. The latter 
is again the well-studied f/(l) s-wave superconductor. Its conductivity has been already 
calculated in [;3j. To check our numerics we have re-calculated it and in Figure 11 we show 
its behavior. It coincides completely with p^. The real part shows the a; = delta function 
characteristic of super conduct ivitjj^ Numerically this can be seen through the l/w behavior 
in the imaginary part. The Kramers-Kronig relation (see (121) in appendix |A]) implies then 
infinite DC conductivity. The real part of the AC conductivity also exhibits a temperature 
dependent gap. 




Figure 11: Real part (left) and imaginary part (right) of the conductivity as a function of 
frequency. The plots correspond to temperatures in the range T/Tc ~ 0.91 — 0.41, from red 
to purple. As expected, the plots reproduce the ones of 



4.5 Conductivities in the (1) — (2) sector 

The relevant equations for the (1) — (2) sector read 



= faT + f'a?^ + 



(2) 



1 



92 

/ 



These equations obey the symmetry 

(a 



(97) 
(98) 

(99) 



®In general, this behavior is also typical of translation invariant charged media, in which accelerated 
charges cannot relax. However, working in the probe limit we effectively break translation invariance and 
therefore the infinite DC conductivity is a genuine sign of superconductivity. 
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One can see that the fact that 0(1) = imphes that a^\l) is independent of a 



.(2)/ 



SO, 



after imposing infaUing boundary conditions at the horizon, the parameter space of boundary 
conditions is two-dimensional, as expected. 

In the unbroken phase the system completely decouples 



(100) 



& ^22 



4.5.1 Diagonal Conductivities a 

The diagonal components of the conductivity, cr^^ and have the same behavior, as 



could be anticipated from the equations (97), (98). Henceforth, we will only refer to a^^, but 
all the conclusions also apply to cr^^. 



Figure 12 shows the conductivity for several values of the temperature. We find that 
these conductivities also show delta-function singularities at w = 0. 




Figure 12: Real (left) and imaginary (right) parts of a^^ versus u for five different temper- 
atures chosen in a range T/Tc ~ 0.91 — 0.41, from red to purple. Im(cr^^) clearly blows up 
as — J- 0. 



The strength of the delta function can also be computed. It is given by the residue of 
the imaginary part of the conductivity at w = 0, 



lim u Im(cr^^) ~ 



(101) 



The residue is plotted in Figure 13 as a function of T/T^. As expected, it starts growing 



from a zero value. At T/T^ ~ 0.65, Ug reaches a maximum and starts decreasing fast, 
changing sign at T/Tc = 0.49. To study Us down to very low temperature we would need to 
go beyond the proble limit. However, as we will comment below, this behaviour of can 
be understood in light of the QNM spectra. 
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Figure 13: Residue at a; = as a function of T/Tc 



Let us look in detail at the behavior of the real part of the conductivity (left plot in figure 
12). For high enough temperatures the optical conductivity is almost constant, Re(cr^^) = 1, 



which is expected since in that regime the dynamics is described essentially by (100). As 



soon as we decrease the temperature, the onset of the DC conductivity also decreases and 
only approaches the constant value asymptotically, when u becomes lar ge en ough and thus 



the term y dominates, turning equations (97), (98) approximately into (100). According to 



the Ferrell- Glover sum rule, the area missing as we lower the temperature is proportional to 

Interestingly, at low temperatures the real part of a^^ starts developing a bump at small 
values of a; (0 < w < 2). The bump leaves less area for the delta function to cover, which 
explains why Ug starts decreasing approximately at this temperature. Moreover, the appear- 
ance of these bumps can be traced back to the fact that for a subleading gauge QNM with 
small I Im(co')|, Reci;(T) >> Imu;(T) holds. Hence, the conductivities affected by this mode 
start developing the reminiscence of a resonance at a particular frequency. We have studied 
the spectrum of low lying QNM for the gauge sector and found that this mode corresponds 
in the normal phase to the lowest excitation of a/i , w = But it is at lower temper- 

atures where one finds a remarkable fact: at T/Tc ~ 0.395 the mode becomes unstable, and 
indeed, as we will see, several physical quantities modify their behavior at that temperature. 

Therefore, we expect a new phase transition around T/T^ ~ 0.395, due entirely to the 
(1) — (2) sector. Since this phase transition seems to be triggered by an unstable mode in 
the vector sector it most likely leads to the formation of a p-wave condensate. This surely 
requires further investigation that is currently underway [55]. 
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4.5.2 Off-diagonal conductivities a & a 

The off-diagonal elements of the conductivity matrix are also related via the symmetry 

-cr^^. Therefore, it is enough to comment on o"^^, although 



(99) and therefore obey a 
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the conclusions are valid for both components. 



The form of is plotted in Figure for various different temperatures as a function 



29 



of frequency. At T/T^ = 1 the system is practically decoupled, so for all temperatures the 
off-diagonal conductivity goes to zero as u increases. 




Figure 14: Real (left) and imaginary (right) part of o"^^ as a function of u for T /T^. ^ 
0.91 — 0.41, from red to purple. 

Observe that a^'^{(jj) behaves as a normal conductivity. Its real part vanishes as w — )■ 0, 
whereas the imaginary part tends to a constant value. 

4.5.3 Conductivities (j+ and a ^ 

It is worth to notice that the equations (97)-(98) decouple if we define a new vector field 



A. 




= S (p 



(102) 



In this basis, the equations of motion become 



= fAl + f'A'^ + [^^^^^^ - ] A^. 



(103) 



It is easy to check that the relation between the conductivity matrices in the two basis is 
given by 

^=^s^y\s-\ (104) 

and that only the off-diagonal components of a are non vanishing. 



The conductivities a ^ and o"+_ are represented in Figure 15 and 16, respectively. The 



plot of the conductivity o"_+ is particularly suggestive. Besides the superconducting delta 
of the DC conductivity, it resembles the behavior observed in Graphene [21]. Such a resem- 
blance of the conductivities of holographic superconductors to the one of graphene has been 
pointed our already in [37] . We emphasize however that the conductivities shown in figure 
15 have an even closer resemblance to [29]. In particular, at small frequencies we see that 
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a Drude-like peak develops. This kind of behavior in metals is usually due to the presence 
of impurities or lattices, whereas in our case, momentum relaxation would be due to the 
non- vanishing expectation value of the charge density operator (Q3) = fiQ. The resemblance 
holds for not too low temperatures. When lowering the temperature, a gap opens up as for 
the (0) — (3) sector. The real part of a+_ shows the same peak already observed for an 
when decreasing the temperature. For temperatures below T ^ 0.49Tc, the pole in the imag- 
inary part of both conductivities changes sign. Of course, it corresponds to the temperature 
at which the residue changes sign. The onset of the DC conductivity at low temperatures 
grows very fast, becoming divergent at T/Tc ~ 0.395. The presence of such a pole in the 
conductivity is related to the appearance of an instability in the spectrum of excitations of 
the gauge field and therefore with a phase transition to another superconducting phase, as 
already discussed. 




Figure 15: Real (left) and imaginary (right) part of the conductivity a ^ for temperatures 

in the range T/Tc ~ 0.91 — 0.41, from red to purple. 



4.6 Quasinormal Modes 

Let us finally study the QNM spectrum in the (1) — (2) sector. This sector contains the 
fluctuations 17, with i = 1,2, therefore in the unbroken phase the spectrum will contain 
two diffusive modes associated with the two gauge fields. The fluctuations of the scalar field 
in the normal phase were already discussed in section |3| Analyzing the quasinormal mode 



spectrum in the broken phase amounts to solving the system of equations (83)-(90). Details 
of the computation can be found in appendix |Bj 
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Re(o-+_) 



10 





Figure 16: Real (left) and imaginary (right) part of the conductivity 
in the range T/T^ ~ 0.91 — 0.41, from red to purple. 



for temperatures 



4.6.1 Type II Goldstone mode 

As expected within the (1) — (2) sector we find a type II Goldstone mode. As in the 
ungauged model for small enough momentum its dispersion relation can be fitted to 



(105) 



Figure [17] shows the dispersion relation for various values of the temperature in the hydro- 
dynamic regime. The quadratic behavior with momentum is apparent. 



Re((j) 
0.012 

0.010 



_,-> T=0.999Tc 
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Figure 17: Plots of Re(a;) (left) and Im(a;) (right) as a function of the momentum. Thick 
lines correspond to data and thin lines to quadratic fit. At T = 0.995 Tc the real quadratic 
parameter B{T) shows a maximum, see Figure 
accuracy. 



18 Relation (105) is fulfilled with high 
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The temperature dependence of B and C is plotted in Figure [TSj The value at T 



Tc is 

given by the same value as in the ungauged model (46) and in fact can also be cross checked 
by calculating the scalar mode dispersion relation in the unbroken phase at T = Tc since 
the QNMs must be continuous through the phase transition. We find a rather surprising 
dependence of B with the temperature. It starts at a finite value at the transition and then it 
rises rather sharply and falls off slower. It reaches a minimum at T ^ 0.49 T^, temperature at 
which we found the change of sign in the residue of current-current correlators. We also find 
another peak around T ^ 0.4 Tc. We expect that it is again related with the instability found 
in the gauge sector around that temperature. It would also be interesting to calculate B(T) 
using an alternative method e.g. as the sound velocity can be calculated from thermodynamic 
considerations alone. In order to do this one would need to formulate the hydrodynamics of 
type II Goldstone modes. We are however not aware of such a hydrodynamic formulation 
and leave this for future research. 
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Figure 18: B (left) and C (right) as a function of T/T^. The zoom-in shows the peak of C 
close to the transition. Furthermore at T ~ 0.4Tc a sharp peak shows up in both coefficients. 
We relate this feature also to the instability arising in the vector sector. 



The attenuation C{T) decreases rapidly with temperature. For temperatures T/Tc < 0.9 
it is negligible and the width of the type II Goldstone scales with in the hydrodynamic 
limit. This fast decreasing with temperature reflects that this mode propagates almost ideally 
in the fluid at low temperature. No further ungapped modes can be found in this sector. 



4.6.2 Higher quasinormal modes 

Higher quasinormal modes correspond to gapped modes in the QNM spectrum and thus 
represent subleading contributions to the low energy Green's functions. We will focus here 
only on two of them: the continuation of the two diffusive modes of the unbroken phase and 
the special gapped mode that appears as the partner mode of the type II Goldstone mode 
in the field theoretical model. 
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Analyzing the first one is interesting in order to understand if also a qualitative change 
in the response pattern, such as that characterized by T* in the U{1) superconductor sector, 
exists in the (1) — (2) sector. Since in this sector there exist however two diffusive modes in 
the unbroken phase it is also possible that the diffusive modes do not simply develop a gap 
but that they pair up and move off the imaginary axis in the broken phase. Indeed as we 
will see this is what happens. 

The special gapped mode corresponds to a mode that is associated to the complex con- 
jugate of the scalar perturbation in the unbroken phase. At = and /i = the scalar 
mode and its complex conjugate are degenerated. As we lower the temperature they split 
into two different modes. When we reach T = Tc, the lowest scalar mode becomes the type 
11 Goldstone mode whereas the mode of the complex conjugate scalar field turns into the 
special gapped mode. The gap of this mode is expected to be given by the tree level result 



(13) m- 



Fate of diffusive modes: As already mentioned, in the (1) — (2) sector we have two de- 
generate diffusive modes in the unbroken phase. When going through the phase transition 
these modes can therefore pair up and move off the imaginary axes such that their quasi- 
normal frequencies develop real parts and lie symmetrically around the imaginary axis. We 
expect therefore that in the low energy limit the dispersion relation takes the form 

uj = T{T)+M{T)k\ (106) 

where both coefficients are complex functions and the second mode is located at u' = —u*. 
Besides, we expect the QNMs to be continuous through the phase transition, which in 
particular means that for T = T^, our pseudo-diffusive modes should match the unbroken 
phase values, i.e. r(Tc) = and A^(Tc) = —i. 



The modes at zero momentum are plotted in Figure 19 We see that indeed the gap 
vanishes as T — )■ Tc, whereas the modes split and develop a real part as we decrease the 
temperature. This last feature is exclusive of the non-Abelian system and thus does not take 
place in the usual U{1) holographic superconductor, where the gap is purely imaginary (see 
[in] and comments above). Close to the phase transition, they present a linear behavior in 
temperature, 

T' 



r(T) = (4.1 - 0.8 i) ( 1 - — ) near . (107) 



The temperature dependence of the coefficient of the momentum in ( 106 ) , (T) , is shown 
in Figure [20| The real part rises very steeply just below the phase transition. The imaginary 
part approaches the unbroken phase value at the critical temperature, i.e. Ai{Tc) = —i, 
as is expected for the pseudo-diffusion modes to continuously connect to the normal dif- 
fusion modes through the phase transition. Notice ImA^(T) decreases when lowering the 
temperature. 

Another check of the fact that the pseudo diffusion modes come from the pairing up of 
the diffusion modes of the normal phase is that their dispersion relation at the phase tran- 
sition matches. Therefore the two diffusive modes are continuous through the transition, as 
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Figure 19: Imco' versus Rew at /c = as a function of the temperature. The shape of the 
figure is compatible with T symmetry, since there are two pseudo-diffusive modes. Having 
Reci;(A; = 0) 7^ is characteristic of the non-Abehan case. 
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Figure 20: Real (left) and imaginary (right) part of Ai{T) as a function of T/T^. As the 
temperature approaches Tc, the value of A^(T) reaches the one prescribed by continuity 
through the phase transition. 
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expected for second order phase transitions, however instead of simply developing an imag- 
inary gap to drop out of the hydrodynamic spectrum as for the usual U{1) superconductor, 
they pair up in two modes that on top of this gap also develop a real part. 

The fact that Re(a;) does not vanish for these modes implies that sufficiently close to 
Tc and in the limit A; = 0, the late-time response of the perturbed state will present an 
oscillatory decay of the perturbations, meaning that, contrary to the U{1) case, there will 
not be a temperature at which the late-time behavior changes qualitatively. 



Special Gapped mode: Seeking for this mode is computationally much more involved. 
Its behavior is characterized by a gap that is proportional to fi. In particular, in [3lj it was 
argued that a type II Goldstone mode is accompanied by a gapped mode obeying co^O) = qfi 
with q being the charge of the corresponding field. In our conventions here we have q = 1- 
So we have to look for a mode with u{k = 0) = /i. Furthermore we expect that it connects 
to the lowest mode of the complex conjugate scalar in the unbroken phase. 



In Figure [21] we depict such mode at zero momentum with respect to the chemical po- 
tential /2 in numerical units. Notice that the mode is continuous at the phase transition, as 
expected. We observe the linear behavior with the chemical potential that is predicted the- 
oretically, at least near ft^- It is very difficult to do the analysis when p, > Q due to the high 
computational power demanded to carry out the computation. The mode shows of course 
also a non- vanishing imaginary part which is due to the dissipation at finite temperature. 
We find that the real part above the phase transition can be approximated by 



Rew = 1.10/i near /ic . 



:io8) 



This result shows a deviation from the conjectured behavior which could nevertheless be 
due to uncertainties in the numerics. Let us emphasize here that the numerics involved in 
tracking this mode through the phase transition were rather challenging. 



Figure 21: Real (left) and imaginary (right) part of the special gapped mode versus the 
chemical potential. We encounter the expected linear behavior with yU. The plot covers both 
the unbroken (dashed line) and the broken (solid line) phases. 
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5 Discussion and Outlook 



The main focus of this work was to estabhsh the existence of type II Goldstone modes in 
the quasinormal mode spectrum of a holographic theory dual to a strongly coupled superfluid 
with U{2) symmetry. 

We studied two models, one in which only the overall U{1) symmetry is gauged in the 
AdS bulk and another in which all the U{2) symmetry is gauged. The most important 
finding is that indeed there exist ungapped excitations represented by quasinormal modes in 
the AdS bulk that show the expected but somewhat unusual quadratic dispersion relation 
of type II Goldstone bosons. 

For the ungauged model this does constitute a surprising result. After all, the field 
theory dual to this model does not contain the necessary conserved currents that would 
correspond to the generators of the global SU{2) symmetry. Standard proofs of the Goldstone 
theorem take the existence of such conserved currents for granted. On the other hand 
it is basically guaranteed that one can construct an effective field theory, a simple Landau- 
Ginzburg type model, that captures the essential dynamics of the light modes, i.e. the lowest 
lying quasinormal modes. Such a model would be essentially given by the field theoretical 
model of section |2] and this guarantees the existence of the type II Goldstone modes. However 
one can expect that such an effective field theory approach can capture only the physics of 
the low lying QNMs but not the higher modes. This is indeed what happens: the partner 
mode of the type II Goldstone mode in the ungauged model does not behave in the supposed 
universal way u = qfi. In contrast the corresponding mode in the gauged model does obey 
this relation approximately and the deviation we found could very well be attributed to 
numerical difficulties and uncertainties that arise in the study of the higher QNMs. 

One rather interesting perspective on the ungauged model opens up if we vary the masses 
of the scalar fields in the AdS bulk. If the masses are slightly different, then at the critical 
temperature only one of the two scalars will feature an ungapped QNM (the one with smaller 
mass). The lowest scalar mode of the second one will still be gapped at that temperature. 
As one goes through the phase transition we do not expect this mode to become massless 
at lower temperatures. Rather it should become a pseudo- Goldstone mode with a gap that 
is proportional to the mass splitting. The appearance of the type II Goldstone mode can 
then be interpreted as the effect of a symmetry enhancement at the point in parameter space 
where the masses of the scalars become degenerate. Since this symmetry is not represented 
by bulk-gauge fields we might call it an accidental symmetry. At this point it is difficult to 
resist the temptation to draw a parallel to the conjectured symmetry enhancement of high 
superconductors. In [38] it was suggested that the phase diagram of high Tc superconductors 
can be captured by a unified model with and enhancement of the 5*0(3) x f/(l) symmetry of 
rotations and electromagnetism to a larger SO (5) symmetry. Since high Tc superconductors 
are d-wave rather than s-wave it remains to be seen how our symmetry enhancement mech- 
anism and the resulting type II Goldstone mode can be combined with holographic models 
of d-wave superfluids such as [5l [6] |^ 

^The appearance of unexpected massless modes related to symmetry enhancement in the context of Bose 
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The second model we studied has bulk gauge fields for all of the U{2) symmetry. There 
are several important differences compared to the ungauged model. The most eye-jumping 
one is that now we can also define and study the full set of conductivities corresponding to 
the U{2) symmetry. Nothing special occurs of course in the unbroken phase, there are simply 
four diagonal conductivities for all the four bulk gauge fields. In the broken phase there are 
however interesting new phenomena. In particular there are now off-diagonal conductivities 
that do not simply vanish. In addition we have found that also the diagonal conductivities in 
the (1) — (2) sector, the one containing the type II Goldstone mode, have delta- function poles 
at zero frequency. In this sense this sector is still superconducting. Moreover, going to a 
decouphng basis for this sector leads to a very suggestive result: the conductivity develops a 
Drude-like peak characteristic of metals on top of the infinite DC conductivity. On the other 
hand Landau's criterion for superfiuidity does not hold in this sector. Recall that this says 
that superfiuidity takes place for fiow velocities v that are smaller than the critical velocity 
Vc where Vc = TamiUi{k) / k for all excitation branches i and over all momenta k [lOl HI]. For 
a type II Goldstone mode the critical fiow velocity is clearly zero. 

A second difference concerns the fate of the diffusive modes. In the unbroken phase there 
are simply four diffusive modes, one for each gauge field in the AdS bulk. In the broken 
phase there is one purely imaginary gapped 'pseudo-diffusive' mode in the (0) — (3) sector, 
i.e. in the sector isomorphic to the U{1) s-wave superfluid. Since there is still one unbroken 
f/(l) symmetry there is also a normal diffusive mode for the preserved U{1) symmetry. In 
the (1) — (2) sector we have however two diffusive modes in the unbroken phase. Going 
through the phase transition these two modes can pair up and move off the imaginary axis, 
becoming a pair of usual gapped quasinormal modes with real and imaginary parts in their 
frequencies. Generically the imaginary part of this gap is smaller (i.e. it lies closer to the 
real axis) then the gap of the purely imaginary mode in the (0) — (3) sector. A large, generic 
perturbation will in its late time response pattern excite both the (0) — (3) and the (1) — (2) 
sector. The late time response of the U{2) invariant order parameter \/[0iP~+lO2]^ will 
therefore be dominated by these paired modes and show an oscillatory behavior in contrast 
to the response pattern of the order parameter in the U{1) case [17^ . 

Another remarkable QNM is the special gapped mode, i.e. the partner mode of the type 
II Goldstone boson. At very high temperatures this mode and the one which at T = 
leads to the sound mode are degenerate. As we lower the temperature the gap of these 
modes becomes different and, for T < Tc, it is expected that Re{u{k = 0)) for the Special 
Gapped mode is proportional to qfi [281 EI]- In particular we find u ~ l.l/i even if g = 1 
in our conventions. Unfortunately with the numerical methods employed in this paper we 
found it very difficult to study this mode and the discrepancy can therefore very well be a 
consequence of insufficient numerical accuracy. It is probably worth the effort to study this 
mode with alternative methods such as the relaxation method developed in [l2l |13] . 

There are several generalizations of the U{2) model that seem interesting and could 
be investigated in the future. A straightforward one would be to analyze the p-wave like 
instability commented in section 4 and to look for a stable background at low temperatures. 

condensates was as well found in [5^ 
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The existence of a perturbation that becomes tachyonic means that the system will suffer a 
second order phase transition into a new phase, probably with the remaining U{1) symmetry 
broken. Such an investigation is currently underway. A similar scenario has been found in 
[m HS], in which a gauged version of the field theoretical sigma model undergoes a phase 
transition driven by an anisotropic vector condensate. 

Another possible generalization would be to analyze the model when the backreaction 
onto the metric is taken into account. This introduces the energy- momentum tensor as an 
operator of the dual field theory and thus we expect the usual sound and shear modes to stem 
from bulk metric fluctuations. Moreover, this would allow us to obtain reliable results even 
at very low temperatures and for instance compute the density of superconducting charge 
densities at zero temperature, as well as B(T = 0). 

We have constructed here a simply model with type II Goldstone bosons using a "bottom- 
up" strategy. It is however also interesting to ask if such models can be realized via "top- 
down" D-brane, string theory or M-theory constructions [Ml HZl HH SH l50] . 

Another possible direction of research involves using the Fluid/Gravity correspondence 
[SI] in order to derive the Hydrodynamic expansion of the current and upon including back- 
reaction also the constitutive relation for the energy-momentum tensor. This will throw light 
on the hydrodynamic behavior of non-relativistic superfluids and in particular should result 
in the formulation of the hydrodynamics of relativistic type II Goldstone modes. Up to our 
knowledge this is not even known to the leading, i.e. zeroth order in derivatives. 

Another direction of investigation concerns the Landau criterion of superfluidity. Accord- 



ing to this criterion, the dispersion relation (105) prevents the system from accommodating 



a superflow. Therefore, even though an and (J22 are superconducting, as soon as a non- 
vanishing supercurrent /superflow is switched on the system should in principle be taken 
out of the superfluid phase. It should be noted however that in a holographic superfluid 
the condensate and its flow are of leading order in a large N expansion and the excitation 
spectrum, the QNMs, are subleading. Therefore it seems not clear if Landau's criterion can 
be applied straightforwardly. It is known however that for the U{1) superfluid there exists a 
critical superflow or a critical supercurrent above which the condensate vanishes [52l |53l [M] • 
It would be very interesting to analyze if this is so (along the lines of [51]) and also to study 
how the whole mechanism takes place. Such an investigation is currently underway [55] . 

Finally it is also interesting to ask the question if holographic models featuring Goldstone 
modes with higher order dispersion relation u = ck"' with n > 2 can be constructed. 

A Matrix valued Kramers-Kronig relation 

The generically matrix-valued spectral function is defined as 

p.,(x) = ([O.(x),O,(0)]), (109) 
where Oi are Hermitian operators. Its behavior under Hermitian conjugation is 

p(x)t = p(-x) = -p(x)*. (110) 
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Correspondingly, the Fourier transform p{k) = J d^x e *'^^p(x) also satisfies a set of identities 

~p{ky = m = -~p{-ky . (m) 

In particular this means that the diagonal components are real and antisymmetric under 
k — )■ —k. One may also be interested in the behavior under u —u. We take now 
k = {oj, q). For theories with rotational invariance the spectral function can depend only on 
q^. Consequently the diagonal components will also be real and odd in u 

Pii{u, q^) = pii{u, q^y = -pii{-u, q^) . (112) 

For the off-diagonal components however, only if one also imposes time reversal or parity 
symmetry can one prove that the off-diagonal entries must be either even or odd functions 
of the frequency. In the present case time reversal symmetry is broken by the presence of 
the chemical potential. Further constraints can however by obtained by supposing that the 
theory is invariant under x — > — x. For an odd number of spatial dimensions we could use the 
parity operators P to take x — )■ — x. In the two spatial dimensions we study in this paper we 
can take P to by a rotation by tt (for an arbitrary even number of spatial dimensions D = 2n 
we could take the angle vr for all the rotations in the i,i + 1-th plane for all i < n). This 
P-operator acts as POi{t,x.)P^^ = aiOi{t, — x) with cTj = ±1. In odd spatial dimensions (Xj 
is the parity of the operator. In even spatial dimension cxj = — 1 if (9j is the component of a 
spatial vector. Hence 

P[pijit, x)] = aiajPij{t, -x) . (113) 
P-invariance implies pij{t,x.) = aiajPij{t, — x), which for the Fourier transform implies that 

pij{uj, q) = -aiajPij{-uj, q)* . (114) 

So the off-diagonal entries are either odd or even functions of u depending on the signs cTj. 
In the case where the fields transform in the same way under the parity operator this means 
that the real (imaginary) part of the off-diagonal components is an odd (even) function of 
the frequency. 



From the spectral function, as defined in (109) we can define two causal propagators 



namely the retarded and advanced Green's functions 



Gn{x)=-iQ{t)p{x), (115) 
GAix) = ie{-t)p{x) , (116) 



where x = (t,x). Using (111), one can prove the following relation among the Fourier 
transforms of these 

GR{k)=GR{-kr = GA{k)K (117) 

From here, we see that the real (imaginary) part, Re^Gn) (lm{Gji)), is even (odd) under 
k — )■ —k. We can compute the Fourier transform of the retarded Green's function, which is 
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given by the convolution of the Fourier transform of the Heaviside step function Q{u) with 
the Fourier transform of the spectral function p{k), 



Using the Fourier transform of the step function 



e{uj - /i)p(/i,q' 



'2^ 



'1181 



u + ie 

and the Sokhatsky-Weierstrass theorem we get 

p{uj', q) du' 



Gr{oo,ci)=V 



0(00, a) 



(119) 



where V denotes the principle value. From the Hermiticity of p{k) we see that we can regard 
(119) as a split of G^{k) into its Hermitian and anti-Hermitian parts, and find that the 
spectral function can be computed from the anti-Hermitian part of the Fourier transform of 
the retarded Green's function 



p{k) = ilGnik) - GR{ky] = 2iG^^\k) 



(120) 



where the (A) stands for anti-Hermitiaij^ Plugging this back into (119) and taking the 
Hermitian part {H) on both sides we arrive at 



G 
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G 



uo — 00' 



-duo' . 



[121] 



which is nothing but the Kramers-Kronig relation for the matrix Green's function. It is 
complemented by the conjugate relation interchanging the Hermitian and anti-Hermitian 
parts. Imposing P-invariance and using (114) and (119) if follows that the Green's function 
satisfies 

Gg(u;,q) = a,a,Gg(-a;,q)*. (122) 

This constrains the QNM spectrum. Taking for example a diagonal Green's function with 
i = j and writing it as a sum over quasinormal frequencies [321 ES] one seems that the 
quasinormal frequencies have to come either in pairs obeying Un and £D„ = — cj* or are 
confined to lie on the imaginary axis. The residues of the pairs are related by complex 
conjugation and the purely imaginary ones have to have also purely imaginary residue. 



B Solving the fluctuation equations 

The (1) — (2) sector of the gauged model in the broken phase consists of a system of 
coupled equations (83)-(90). In order to extract the spectrum of quasinormal modes we 



Using (117) we can always work with retarded Green's functions Gr. 
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made use of the techniques detailed in [TSl El], where a method to compute the poles of 
the Green functions in terms of non-gauge invariant fields was developed. The quasinormal 
frequencies are given by the zeroes of the determinant of the field matrix spanned by a 
maximal set of linearly independent solutions satisfying infalling boundary conditions on the 
horizon evaluated at the boundary. 

Imposing infalling boundary conditions, the near horizon behavior of the fields solving 
the mentioned equations reads 



a = (p-l)'^(«(o) + «(i)(p-l) + ...) , (123) 

/3 = (p-l)-(/3(o) + /3(i)(p-l) + ...) , (124) 

a« = (p - 1)-+^ (a%) + a« )(p - 1) + . . . ) , (125) 

0^ = {p-ir[a%^ + af^,^{p-l) + ...) , (126) 



where k = —iu/3 and i = 1,2. Since the system is subject to two constraints, we can only 
choose four of the six parameters at the horizon. Without loss of generality, solutions can be 
parametrized by {a(o), /3(o), a^*(o)}- ■'■^ ^^^^ possible to construct four independent 

solutions to the field equations. We can label them as /, //, III, IV. 

Two additional solutions, V, VI, can be obtained by performing gauge transformations 
of the trivial solution, 

13^ i^, ^ -kXi, af'^ 0, aj^^ -> uXi, af ^ ^ iQXi , (127) 
« ^ /9 ^ 0, ^ 0, ^ -kX2, af ^ ^ -iGAs, af ^ ^ cjAa , (128) 

where Aj are arbitrary constants. Notice that these pure gauge solutions are not algebraic 
since they have a nontrivial dependence on the bulk coordinate p. 

The most general solution for each field yjj = {a, /3, cl]^ , ai*''} is given by a linear combi- 
nation of the above solutions, including the pure gauge modes, 

Lfi = Ciifl + Cn^l^ + CiiM" + C/yV^i^ + Cy^jf + Cy/^jf ^ , (129) 

where we have defined {a;(p), /3(p)} = {pa(p), p/3(p)}. This convenient choice allows us to 
identify the asymptotic boundary values ipi with the sources of the gauge invariant operators 
of the dual field theory. 

As shown in [15|, the poles of the retarded Green functions will be given by the values of 
the frequency for which the determinant of the matrix spanned by iff vanishes asymptoti- 
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cally. Expanding the determinant and evaluating it at a cutoff p = A, it reads 



A1A2 



det 



II 
II 



III 



II 



^t{\) ^t{\) ^t{\) 



III 
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III 



II 
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VP 



VI 



\ 



VI 



IV 
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(130) 



jv 



Va 

v7 

Vil) / 



vi"" 

,jv 



where the background boundary conditions 0(A) = and A\l/ = have been aheady im- 
posed. This absence of background sources for the corresponding operators makes the point 
(cj, k) = (0, 0) a trivial solution to the vanishing determinant condition, which ensures the 
existence of a hydrodynamic mode. Notice also that the point {u, k) = (0, 0) is a double 
solution to the previous determinant equation. 

Solutions to the equations of motion and to the determinant condition (130) have been 
computed numerically. It has been checked that the election of solution basis, i.e. of initial 
values of the free parameters, does not affect the result. 
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